Hadronic effects on the X(3872) meson abundance in heavy ion colHsions 
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We study the hadronic effects on the X(3872) meson abundance in heavy ion collisions. We 
evaluate the absorption cross sections of the X(3872) meson by pious and rho mesons in the hadronic 
stage of heavy ion collisions, and further investigate the variation in the X(3872) meson abundance 
during the expansion of the hadronic matter for its two possible quantum number states; = 1^ 
and 2~. We show that the absorption cross sections and the time evolution of the X(3872) meson 
abundance are strongly dependent on the structure and quantum number of the X(3872) meson. We 
thus suggest that studying the abundance of the X(3872) meson in relativistic heavy ion collisions 
provides a chance to infer its quantum immber as well as its structure. 

PACS numbers: 14.40.Rt, 25.75.-q, 13.75.Lb 



I. INTRODUCTION 

Relativistic heavy ion collision experiments have made 
it possible to study a system of quantum chromodynamic 
(QCD) matter at very high temperature and density in 
the laboratory. The research on the system of decon- 
fined quarks and gluons, so called quark-gluon plasma 
(QGP), has also enabled us to understand the possible 
phase transition between the hot and dense matter and 
QGP. Moreover, due to the enormous energies produced 
in heavy ion collisions, particles that are otherwise hard 
to find in nature could be produced during the quark- 
hadron phase transition. 

Recently STAR Collaboration reported the observa- 
tion of an antimatter helium-4 nucleus as well as an an- 
timatter hypernucleus produced at RHIC [11, [^j , and also 
tried to measure the signal of an exotic H dibaryon Q. 
As attempts to understand the production of particles 
of these kinds, there have been many studies focusing 
on their production yields based on both the statistical 
model and the coalescence model Q-Q- Moreover, as 
one of the possible methods to understand the structure 
of the exotic hadrons, a new approach of studying exotic 
hadrons in relativistic heavy ion collision experiments has 
been proposed 0,01 • There, the relation between the pro- 
duction yields of exotic hadrons and the structure at the 
moment of their formation has been sought out by consid- 
ering the production of all proposed possible structures 
using the coalescence model, and it was found that the 
production yields of exotic hadron candidates strongly 
reflect their structures. 

The abundance of hadrons evaluated at the chemical 
freeze-out temperature, however, may change due to the 
dissociation or the absorption by mostly light mesons 
such as the pion and the p meson in the hadronic medium. 
The effects from the hadronic interactions on the pro- 
duction of heavy quark mesons have been discussed in 
many literatures. In order to estimate the possibilities 
of J/tP suppression in the hadronic matter, one meson 
exchange model with the effective Lagrangian has been 
introduced to evaluate the absorption cross sections with 
light hadrons [lol - [l3j . A similar approach has been ap- 



plied to investigate the time evolution of Dsj(2317) me- 
son abundance in the hadronic matter [l^ . In this work 
we investigate the hadronic medium effects on the pro- 
duction yield of one of exotic mesons, the X(3872) meson. 

The X(3872) meson was first discovered by Belle 
Collaboration [l5| from the measurement of — 
J/\['7r+7r~/C+ decay, and later confirmed by CDF jT6| . 
DO [13, and BABAR IHl CoUaborations. The addi- 
tional decay modes of X(3872) mesons to D°D°7r° [H- 
[U, J/^uj [13, J/^-T [ll, and *(2s)7 [H have also 
been observed. The positive charge parity of the X(3872) 
meson has been established by the observation of the 
X(3872) meson decaying to J/^j [H and 1'(2s)7 [H, 
and the current world average mass of the X(3872) meson 
in PDG [H is 3871.68 ± 0.17 MeV. However, it is still 
not clear what the exact structure and quantum number 
of the X(3872) meson is. Suggested hypotheses for the 
structure of the X(3872) meson include a pure charmo- 
nium state, a D^D*^ hadronic molecule, a tetra-quark 
state, and a charmoniun-gluon hybrid state [26| . From 
the analysis of the angular distribution of the X(3872) 
meson decaying to J/^'tt+tt" [27j . we now understand 
that the possible quantum number should be either 
1+ or 2". 

There are various experimental results supporting each 
spin possibility of the X(3872) meson. The study of the 
X(3872) meson decaying to disfavors the 2~ quan- 

tum number because of the angular momentum barrier 
in its near-threshold decay [l^H^. On the other hand, 
the analysis of the X(3872) meson decaying to J/vtw fa- 
vors a p-wave in the final state, 2~ [13 . We expect that 
the two different spin possibilities of the X(3872) meson 
will also lead to different experimental results in heavy 
ion collision experiment. 

After the X(3872) meson is produced at the chem- 
ical freeze-out, it interacts with other hadrons during 
the expansion of the hadronic matter. As a result, the 
X(3872) meson can be absorbed by the co-moving light 
mesons or additionally produced from interactions be- 
tween charmed mesons such as D and D* . Evaluating 
the X(3872) meson cross sections by light hadrons there- 
fore should be useful in estimating the hadronic effects 
on the X(3872) meson abundance in heavy ion collisions. 
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However, the X(3872) meson would interact with hght 
hadrons differently depending on the spin of the X(3872) 
meson. In order for the spin-2 X(3872) meson to inter- 
act with light mesons, there should be a charmed me- 
son having the relative momentum to satisfy the angular 
momentum conservation. Also the spin-2 X(3872) me- 
son should carry a symmetric traceless spin polarization 
tensor whereas the spin-1 X(3872) meson carries a po- 
larization vector. Therefore, we expect to obtain two 
different results when we evaluate the cross sections of 
the X(3872) meson for the two different spin states. By 
comparing these results with the experimental observa- 
tion in heavy ion collisions, we may obtain a hint for the 
quantum number of the X(3872) meson. 

In this study we restrict our consideration of the 
X(3872) meson structure to the spin-I tetra-quark state 
and the spin-2 pure charmonium state. We briefly dis- 
cuss the D^D*° hadronic molecule for the spin-1 X(3872) 
meson possibly produced during the hadronic stage and 
at the kinematical freeze out point. All the discussion 
will be focused on the central heavy ion collisions at Rel- 
ativistic Heavy Ion Collider (RHIC) at Brookhaven Na- 
tional Laboratory; using a model developed to describe 
the dynamics of the cental Au-|-Au collisions at ^/snn = 
200 GeV. Hereafter, we use simplified notations for the 
X(3872) meson; Xi for a 1"*" state and X2 for a 2~ state. 

This paper is organized as follows. In Sec. II, we 
briefly discuss the production of the X(3872) meson at 
the chemical freeze-out in both the statistical and coa- 
lescence model. Then we consider the interaction of the 
X(3872) meson with light mesons such as pions and p 
mesons, and evaluate the cross sections of the X(3872) 
meson in the hadronic medium in Sec. HI. In Sec. IV 
we investigate the time evolution of the X(3872) meson 
abundance by solving the kinetic equation based on the 
phenomenological model. Sec. V is devoted to conclu- 
sions. In Appendix A, we discuss the dependence of the 
strong coupling constants on the X(3872) meson mass. 



We consider the X(3872) meson produced at the 
hadronization temperature = 175 MeV when the vol- 
ume of the quark-gluon plasma, Vh is 1908 fm'^. We 
assume that the total number of charm quark produced 
from the initial hard collisions at RHIC is 3, which leads 
to the charm quark fugacity factor 7c = 6.4 by the 
requirement that the charm quark is conserved among 
charmed hadrons such as I?, D*, Dg mesons, and Ac. 
The difference in the yields between the spin-1 X(3872) 
meson and the spin-2 X(3872) meson originates only from 
the spin degeneracy gx in the statistical model. 

In the coalescence model, which successfully explains 
both the enhancement of the baryon to meson ratio in the 
intermediate transverse momentum region j32H35i and 
the quark number scaling of the elliptic flows [36|, we 
consider the yields of the X(3872) meson produced from 
both the four quark configuration for the spin-1 state and 
the two quark configuration for the spin-2 state. We as- 
sume that the quark coalescence occurs in the volume 
1000 fm'^ and the mass of light constituent quarks is 300 
MeV, while that of a charm constituent quark is 1500 
MeV. We also assume that the available light quark num- 
ber at hadronization temperature is 245. We adopt the 
oscillator frequency of the Wigner function for charmed 
hadrons = 385 MeV obtained by requiring that the 
coalescence model reproduces well the yield of Acf2286) 
in the statistical model. For details, refer to Ref. [5[- We 
summarize the production yields of the X(3872) meson 
in Table n 



TABLE I: The X(3872) meson yields in central Au-|-Au col- 
lisions at y/SNN = 200 GeV at RHIC in both the statistical 
and coalescence model. The yields for the spin-2 state differ 
from those in the Table V. of Ref. by the spin degeneracy 
factor 5/3. 



II. X(3872) MESON PRODUCTION FROM THE 
QUARK-GLUON PLASMA 

We evaluate the production yields of the X(3872) me- 
son in heavy ion collisions using both the statistical and 
the coalescence model. The statistical model, which as- 
sumes that hadrons are in thermal and chemical equilib- 
rium when they are produced at chemical freeze-out in 
heavy ion collisions, has been very successful in describ- 



ing the production yields of hadrons [2S 
the same parameters evaluated in Ref. 
thermal yields. 



M. We apply 
5| to obtain the 
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In Table U the smaller yields in the coalescence model 
compared to those in the statistical model reflect the sup- 
pression effects in the quark coalescence process. The 
suppression mechanism is, however, different for differ- 
ent spin states of the X(3872) meson. The coalescence 
of additional quarks to construct the tetraquark state 
makes the yield suppressed for the spin-1 X(3872) meson 
whereas the construction of a d-wave coalescence factor 
leads to the smaller yield for the spin-2 X(3872) meson. 
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FIG. 1: Born diagrams for X(3872) absorption by pions and p mesons. X-k — >■ D*D*; (a) and (b), Xtt — >■ DD; (c) and (d), 
Xp^ D*D; (e) and (f), Xp-t 0*0- (g) and (h), Xp DD; (i) and Xp^ D*D'; (j). 



III. HADRONIC EFFECTS ON THE X(3872) 
MESON 

The X(3872) meson produced at the chemical freeze- 
out interacts witli other hadrons during the expansion of 
the hadronic matter. As a resuh the X(3872) meson can 
be absorbed by the co-moving hght mesons or produced 
from the interaction between charmed mesons such as D 
and D*. We consider here the X(3872) meson interacting 
with hght mesons such as pions and p mesons. 
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D*D, 
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The diagrams representing each process in Eq. ^ are 
shown in Fig. [TJ To evaluate the cross sections for these 
diagrams in Fig. [1] we consider the following interaction 
Lagrangians 



C-^DD' igTrDD'D*^'T ■ {Dd^n - d^Dn) +II.C. 

CpDD = igpDDiDfdf^D - dpDrD) ■ p^ , 

CpD'D' = tgpD'D' [{d^D*''TDt-D*''Td^Dt) 

+ {D*''t ■ d^p, - d^D*-f ■ p,)D~*t^ 

+ D*^^iT-ffd,Dt-r-d^frDt)], 

C^DD = ig^DDi^^{Dd^D - dy,DD), 

C^D'D' = ig4,D'D'[V{d^.D*''Dl-D*''d^Dl) 

+ id^rD:-vd^D:)D~*^ 

+ D*'^ird^Dt~d^rDt)]. 



(3) 



In Eq. T arc the Pauli matrices, and if and p de- 

note the pion and rho meson isospin triplets, respectively, 
while D = {D^,D+) and D* = {D*",D*+) denote the 
pseudoscalar and vector charm meson doublets, respec- 
tively. Here the shorthand notation ijj has been used for 
the J/iJj meson. These interaction Lagrangians have been 
obtained from the free Lagrangians for pseudoscalar and 



vector mesons by introducing the minimal substitution 
p^ . On the other hand, the interaction Lagrangians 
for the X(3872) meson have been built to produce strong 
transition matrix elements for the X(3872) meson decays; 
X-^ J/iPp and X^ D^D*° [13 • 



C,XiD*D = gXiD*DXiD*^D, 

(^Xi-ipp = igXitpp^'^'^'^'^i'uPpdcrXlfj, 

t^X2D*D = —igx2D'DX2'^ D*^d^D, 

Cx2-il>p = ~gX24'p<^'"'^'^^2^ia{du'fp°'dpPa - d^p^'dpiJa) 

+ g'x^^pe^^'^d^X^^a^id^^pP.-^Ppd^p,). (4) 

As we see in Eq. we need an additional derivative 
for the interaction Lagrantians of the spin-2 X(3872) me- 
son compared to those of the spin-1 X(3872) meson to 
satisfy the angular momentum conservation. The struc- 
ture of these Lagrangians is expected to make the en- 
ergy dependence on the cross sections different in the 
hadronic medium. This is also the factor prohibiting the 
spin-2 X(3872) meson from decaying to the vector meson 
D* and the pseudoscalar meson D near the threshold en- 
ergy. We need an anti-symmetric tensor e^'^^'^ to describe 
the isospin violating interaction between one axial vector 
meson and two vector mesons for the spin-1 X(3872) me- 
son. This term becomes more complicated for the spin- 
2 X(3872) meson since the anti-symmetric tensor e^'^'"^ 
should not be fully contracted with the symmetric polar- 
ization tensor tt^j^. 

Based on the above effective Lagrangians. we consider 
the reactions for the X(3872) meson absorption by pions 
and p mesons shown in Fig. [1] Among those diagrams 
in Fig. [l]the process Xp ^ D*D leads to the same cross 
section as the process Xp — >■ D*D. The amplitudes for 
all processes, without isospin factors and before summing 
and averaging over external spins, are represented by 
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where the amphtudes for the first process Xir — > D*D* 



My^ = -gTTD'DgXiD'D<^i<^2'^ 
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(6) 



for the spin-1 X(3872) meson and 



^xi ^ -9^D' 09X20*0-^1 £2^3 



f 



-{P3-P2)a{'2p2-P3)i^- (7) 



for the spin-2 X{3872) meson. Similarly, the amplitudes 
for the second process Xtt — > DD are 
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and 
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Finally the amplitudes for the processes Xp ^ DD and 
Xp D*D* are 

^Xi = 9-4>oo9Xi,pp£^''''''eipe2p—^-^{P4-P3)aPia 



(Pl +P2)i/(Pl +P2)" 



■^X2 = ~9i'Oo9x2->pps'^'''''' -^1^0,^2. 
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(12) 



and 



for the 1+ state and the 2 state, respectively. And the 
amplitudes for the process Xp D*D are 
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(13) 



respectively, for both the Xi(3872) meson and the 
X2(3872) meson. 

In the above equations, pj denotes the momentum of 
particle j. We choose the convention that particles 1 and 
2 represent initial-state mesons, while particles 3 and 4 
represent final-state mesons on the left and right sides of 
the diagrams, respectively. Here we use the usual Man- 
delstam variables given by s = (pi -f ^2)^, i = (pi — Pa)^, 
and u = (pi — p4)^. The polarization tensor tt^'^ satisfies 
the following polarization sum. 



incoming particles for s channel taken in the center of 
mass frame. For the cutoff parameter A, we use A = 2.0 
GcV. The final isospin- and spin-averaged cross section 
is given by, 



\Pf\ 



647r2sgi52 \pi 



dn\M\'^F'' 



(16) 



with gi and 92 being the degeneracy factors of the initial 
1 and 2 particles, (2/i-hl)(2S'i + l) and (2/2 + l)(2S'2 + l), 
respectively. We denote by |A^p the squared amplitude 
of all processes in Eq. ([5]) obtained by summing over the 
isospins and spins of both the initial and final particles 
after killing all unphysical terms satisfying pi^ir^'^ = 0, 

P2p£2 = ^'^'^ so on. In evaluating |A^p we only con- 
sider the X(3872) meson interacting with I?" and D^* 
(or and ) since the X(3872) meson does not have 
its isospin partner. In Eq. (|16p . \pi\ and \pf \ represent 
the three-momenta of the initial and final particles in 
the center of mass frame. For coupling constants, we 
use gpDD = gpD*D'=2.52, g^oD = g^D*D*^T-Q^ from 



121. and g-no* D—^-'i from the decay width of D* meson 
38|. The strong coupling constants for the X(3872) me- 
son have been taken from the Table II in Rcf. [33], and 
those are summarized in Table |lll 



TABLE II: The strong coupling constants for the X(3872) 
meson 13711. 



pol 

1/ uu' k^'k'''\( a'v fc^'^'^A 
2 V mjj- / V mj^. / 

= \ [g^^'g"""' + g'-'^'g^''' - g^^'g^'"') - ^ 

X [g^^'K'k''' + g'^'^'k^k^'' -t- g^'^'k^'k'' ^ g^'^'k^k'''^ 

+ \ + ^^''fc') W"' + ^fc^'fc"') • (14) 
6 V m^j. / V mjj- / 

In obtaining the full amplitudes we introduce the fol- 
lowing form factors at interaction vertices to prevent the 
artificial growth of the tree-level amplitudes with the en- 
ergy, 

where is the squared three-momentum transfer for t 
and u channels and the squared three-momentum of the 



y = 1+ P = 2- 



gXjD'D 


3.5 ± 0.7 GeV 


189 ± 36 


gxji,p 


0.14 ± 0.03 


-0.29 ± 0.08 GcV"' 


9xji,p 




0.28 ± 0.09 GoV"' 



In Fig. [21 we show the cross sections for the absorption 
of both a Xi(3872) meson and a X2(3872) meson by pi- 
ous and p mesons via processes Xtt — > D*D*, Xtt DD, 
Xp D*D, Xp DD, and Xp D*D* as functions 
of the total center-of-mass energy s'/^ above the thresh- 

1 /2 I 

old energy Sq of each process. We see in Fig. [2] (a) 
that there exists a peak near the threshold energy for 
the endothcrmic process Xitt — )■ D*D*, while the cross 
sections for the other exothermic processes Xitt DD, 
Xip D*D, Xip -)> DD, and Xip D*D* become 
infinite near the threshold. The same general behaviors 
as these can be found for the spin-2 X(3872) meson in 
Fig. [21(b). For the endothcrmic process X27r — > D*D*, 
however, there is a gradual decrease after a strong rise 
near the threshold. It is also noticeable that there ex- 
ists a sharp dip very near the threshold and a gradual 
increase and decrease afterwards for the exothermic pro- 
cesses X2P — >■ DD and X2P — >■ D*D* , which is similar to 
that shown in Ref. [131 for the dissociation cross section 
of J /tp into open-charm meson mediated by the spin-2 
X(3872) meson. These differences are due to two dif- 
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FIG. 2: Cross sections for the absorption of (a) a Jfi(3872) 
meson and (b) a X2(3872) meson by n and p mesons via 
processes Xvr D*D*, Xn DD, Xp D'D, Xp DD, 
and Xp D*D*. 



ferent interaction mechanisms originated from two possi- 
ble spin quantum numbers. The additional derivative in 
the interaction Lagrangian for the spin-2 X(3872) meson 
causes the completely different energy dependence on the 
cross sections especially for the processes Xp — DD and 
Xp^ D*D*. 

We also clearly see in Fig. [2] that the absorption cross 
sections for the processes Xtt — > D*D*, Xtt — >■ DD, and 
Xp — )• D*D are much bigger when the spin of X(3872) 
mesons is 2 than when it is 1. This is largely attributed 
to the large strong coupling constant used to evaluate the 
cross sections for the processes having X2D*D interac- 
tions. For the processes involving the p meson, the ad- 
ditional derivative in the interaction Lagrangian for the 
spin-2 X(3872) meson causes 10 to 40 times bigger cross 
sections in the processes X2P — > DD and X2P — >■ D*D* 
compared to those in the the processes Xip — >■ DD and 
Xip — J- D*D* as shown in Fig. [2] when the strong cou- 



plings gxi-4,p, gx^ipp, and ^3-^^ ^^^'^ ^^^'^^ 

processes Xtt D*D*, Xn' DD, and Xp D*D, 
however, the additional derivative brings out roughly 
Tod ^ 1.9 GeV. This factor will multiply the already 
large coupling constant gx2D'D given in Table HIl mak- 
ing the effective coupling strength ( ^ 189 x 1.9 GeV ) 
much larger compared to gxiD'D = 3.5 GeV. This ex- 
plains the bigger cross sections for X2(3872) meson than 
those for the Xi(3872) meson in Fig. [1 

Both strong coupling constants gxiD*D and gx^D'D 
were obtained from one experimental measurement us- 
ing two different spin possibilities [13] ■ As was already 
pointed out, however, the analysis of the X(3872) meson 
decaying to disfavors the 2^ quantum number 

because of the angular momentum barrier in its near- 
threshold decay [l^, [l^] ■ The D meson should have the 
relative angular momentum in order to be able to interact 
with the spin-2 X(3872) meson to satisfy the angular mo- 
mentum conservation. Therefore, the coupling constant 
gx2D*D has to be large to compensate for the angular 
momentum suppression near threshold. 

Nevertheless, it is still possible to get a smaller strong 
coupling constant gx2D' d when the X(3872) meson mass 
increases slightly. In Appendix A, we have investigated 
the origin of the big strong coupling constant gx2D'D^ 
and have found that it is very sensitive to the variation 
of the X(3872) meson mass. Varying the mass within 
the experimental uncertainty, we estimate that 5x2-0*1) 
could be reduced by a factor of -\/3 and the cross sections 
of the spin-2 X(3872) meson evaluated in Fig. [5] by a 
factor of 3. 

The bigger cross sections for the spin-2 X(3872) meson 
are contrary to naive expectations. It is expected that 
the size of the bag containing four quarks should be at 
least bigger than that of the bag having two quarks. In 
the simple bag model the size of the bag increases with 
the number of quarks inside the bag as i? oc N^^^ [39} . 
Since the cross section depends on the size of the hadron 
in general, we expect the cross section of the X(3872) 
meson composed of four quarks to be bigger than that 
of the X(3872) meson made up of two quarks. However, 
we sec here only the effects from the interaction mecha- 
nism caused by two different spins since the interaction 
Lagrangians arc blind to the size of the X(3872) meson. 



IV. TIME EVOLUTION OF THE X(3872) 
MESON ABUNDANCE IN HADRONIC MATTER 

Using the cross sections evaluated in the previous sec- 
tion we now consider the time evolution of the X(3872) 
meson abundance in hadronic matter. We build the 
evolution equation consisting of the densities and abun- 
dances for hadrons participating in all processes shown 
in Fig. [T] TT, p, D* , and D mesons. 
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FIG. 3: Thermally averaged cross sections for the absorption of (a) a Jri(3872) meson and (b) a X2(3872) meson by pions 
and p mesons via processes Xtt — ^ D* D* , Xtt — >■ DD, Xp — > D*D, X/o — ^ 7? I), and Xp — >■ D* D* and their inverse processes 
->X7r, I)D ^Xtt, D'D ->Xp, DI) ^Xp, and D' D* ^Xp for (c) a Xi(3872) meson and (d) a X2(3872) meson. 



- i?QGp(r) + [{<^cc'^IXVcc'W{t)NAt) 
-{crix^cc'Vix)ni{T)Nx{T)^, (17) 

where ^/(t) and nc(T) are, respectively, the density of a 
light meson such as a pion or a p meson and the density of 
a charmed meson in the hadronic matter at proper time 
T, whereas Nc'{t) is the abundance of the other charmed 
meson in each process shown in Fig. [T] at proper time 
T. ni{T), ndr) and Nc'{t) are calculated from Eq. ([T]) 
by assuming that light mesons and charmed mesons are 
in equilibrium and vary in time through the temperature 
profile introduced below, Eq. p9)) . In the above rate 
equation, Eq. ([T7|) . {(Jab^cdVab) is the cross section aver- 
aged over the thermal distribution for initial two particles 
in a two-body process ah — )■ cd given by plo| 



{c'ab^cdVab) 

_ J d^Pad^Pbfa{Pa)fb{Pb) O" ab^cd'Vab 
J d^Pad^Pbfa{Pa)fb{Pb) 
1 

= A 9r- r ^ 2T- r ^ / dzKi{z)a{s = z^T^) 
^alK2{aa)alK2{ab) A„ 

X [z^ - [ua + abf] [z^ - {aa -abf], (18) 

with Oii = rrii/T^ zq = max(aa -|- ab,ac + ctd), Ki and 
K2 being the modified Bcssel function of the first and 
second kind, respectively and Vab denoting the relative 
velocity of initial two interacting particles a and 6, Vab = 
\/iPa -PbY - mlml/{EaEb). 

The X(3872) meson abundance at proper time r, 
Nx{t), depends on both the dissociation rate such as 
Xtt -> D*D*, Xtt DD, Xp ^ D*D, Xp DD, 
and Xp — )• D*D* and the production rate through the 
inverse processes, D*D* — j-Xtt, DD — j-Xtt, D*D ^-Xp, 
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DD —>-Xp, and D*D* — >-Xp. We use the detailed balance 
relations based on the results for the forward processes 
shown in Fig. [2] in evaluating the thermally averaged 
cross sections of the inverse processes. The results are 
shown in Fig. [3] 

In Eq. (|17|) . 71;(t), ndr) and Nc'{t) varies in time 
through the temperature profile developed to describe 
the dynamics of rclativistic heavy ion collisions. We 
use the schematic model based on the boost invariant 
Bjorken picture with an accelerated transverse expansion 
[ij, Hll . The system of the quark-gluon plasma of its fi- 
nal transverse size Rc at the chemical freeze-out time tc 
expands with its transverse velocity vc and transverse 
acceleration ac- The temperature of the system is main- 
tained with a constant temperature Tc until the end of 
the mixed phase at th, and decreases afterwards to the 
kinetic freeze-out temperature Tp. The volume and tem- 
perature profiles as a function of the proper time r are 
as follows. 



V{t) = tt[Rc + vc{t - Tc) + ac/'2{T - tcY 



T{t) ^Tc^ (Th - Tp) 



Tp - Th 



4/5 



(19) 



with Th and Tp being the hadronization temperature, 
and the freeze-out time, respectively. The values used in 
Eq. ([T^ are summarized in Table Hill 



TABLE III: Values for the volume and temperature profiles 
in the schematic model Eq. (|19p . 



Temp.(MeV) Time (fm/c) 



Rc = 8.0 fm Tc = 175 
vc = 0.4c Th = 175 
ac = 0.02cVfm Tp = 125 



TC = 5.0 
T„ = 7.5 
TP = 17.3 



In solving Eq. ([T7| we have assumed that the total 
number of charm quarks in charmed hadrons is conserved 
during the evolution of the hadronic matter. It has been 
discussed that chances for charmed mesons to be pro- 
duced and destroyed in the hadronic matter are very 
small because of their small production and annihilation 
cross sections [H], |TJ |42| . The light mesons are assumed 
to be in equilibrium with the medium and the total num- 
ber of the pion is set to 926 at freeze-out [4l| and that 
of the p meson to 68 after considering the contributions 
from the decays of resonances. To take into account the 
effect of the production of the X(3872) meson through 
hadronization from the quark-gluon plasma, we include 
the term Rqgp{t) [13] given by 



TC), Tc <T <Th 

otherwise 



with iV^ being the total number of the X(3872) meson 
in Table |T] produced from quark-gluon plasma either by 
the two-quark coalescence or by the four-quark coales- 
cence. We assume here that the volume of the quark- 
gluon plasma decreases linearly during the phase transi- 
tion time Th — Tc = 2.5 fm/c whereas that of the hadron 
gas increases with a rate enough to occupy both the de- 
creased volume of the quark-gluon plasma and the newly 
increased volume of the entire system by the expansion. 




(20) 



8 10 12 14 16 18 
T(fm/c) 



FIG. 4: Time evolution of X(3872) meson abundances in cen- 
tral Au-Au collisions at ^/sNN ~ 200 GeV for diflFerent states 
produced from the quark-gluon plasma. 

In Fig. m we show the abundances of the X(3872) me- 
son as a function of the proper time for different states 
produced from the quark-gluon plasma in central Au-Au 
collisions at ^/sNN = 200 GeV. Since the scattering cross 
sections for the 1"*" state Xi(3872) meson are so small as 
shown in Fig. [51 the abundance obtained by the four- 
quark coalescence increases very slightly to 4.3xl0~^, 
while the expectation from the statistical model decreases 
also very slightly to 2.8x 10^'*. However, due to the large 
scattering cross sections for the 2~ state X(3872) me- 
son, the X(3872) meson in the normal cc state with d- 
wave has more chances to interact with light mesons in 
the hadronic evolution, and therefore the abundance de- 
creases fast to 1.2x10"^. The thermal model expectation 
for the spin-2 X(3872) meson also decreases rapidly, fol- 
lows the evolution of the coalescence model abundance, 
and cvolutes afterward together. The final ratio of the 
abundance for the X2(3872) meson over that for the 
Xi(3872) meson both in the coalescence model is ex- 
pected to be ^ 2.8 at the kinetic freeze-out. 

Based on the above investigation about the time evo- 
lution of the AT (3872) meson, we can further consider the 
possibility of producing a hadronic molecular state of the 
spin-1 Ar(3872) meson. If the state is a hadronic molecule 
composed of £)°l)*°(Z)°£)*°) in s-wave, it will be domi- 
nantly produced at the end of the hadronic phase through 
hadronic coalescence. The production in the hadronic 



9 



phase through the two body hadronic interaction would 
be very small. For such process to be possible, the in- 
verse processes like p*D* — >X7r, DD -^Xtt, D*D -^Xp, 
DD — J'Xp, and D*D* — >X/9 should prevail the forward 
processes. We have already seen, however, in Fig. m 
that the thermally averaged cross sections of the inverse 
processes are smaller than those of the forward processes. 
However, this does not mean that the X(3872) meson can 
not be produced in the hadronic phase. In fact, hadronic 
coalescence will continue to occur but the total absorp- 
tion should be very large since the size of the loosely 
bound hadronic molecule is thought to be much bigger 
than that of the compact tetraquark state with the same 
spin. 



8.0 



6.0 



4.0 



2.0 



0.0 



-.-d'"''*V (1') 

q^VcC^) 



Stat. (1 



10 12 
T(fm/c) 



14 16 18 



FIG. 5: The possible time evolution of the hadronic molecular 
state of the X(3872) meson produced during the hadronic 
stage in central An- An collisions at ^sjvjv = 200 GeV when 
X(3872) is assumed to be a 1+ state. 

With this in mind we can estimate the number of 
X(3872) mesons in the hadronic phase by solving the rate 
equation, Eq. ([T7| for the hadronic molecular state back- 
wards in time from the final yield of 7.8 x IQ-^ 0, H 
calculated through the hadron coalescence at the end 
of the hadronic stage. As shown in Fig. [Sj when the 
hadronic molecular state of the Xi(3872) meson is pro- 
duced sometime during the hadronic stage by the hadron 
coalescence, the number of X(3872) mesons is expected to 
be in the range between 7.8x10-"^ and 8.1x10"''. Owing 
to the small cross sections evaluated in this work, which 
are blind to the size of the hadron, the yield decreases 
slightly during the hadronic stage, finally resulting to the 
ratio '--^ 18 between the hadronic molecular state and the 
tetraquark state at the kinetic freeze-out. 

The discussion on the production yield for a molecu- 
lar state of the X(3872) meson is based on the assump- 
tion that it is a D'^D*^ state. However, it has been 
shown that charged components of D and D* mesons 
also play an essential role in explaining the branch- 
ing ratio of the X(3872) meson decaying to lo and p 



mesons [43l - l45| . If we take into account a linear com- 
bination of D and D* mesons for the X(3872) meson, 
|X(3872)) = 1/V2(|D0Z)*0) + \D+D*~)), then we have 
to evaluate the average production yields coming from 
both D°D*° and D+D*~. Thus, the production yield 
would be the same since the numbers of D or D* mesons 
are independent of their charges. 



V. CONCLUSION 

We have studied the hadronic effects on the X(3872) 
meson abundance in heavy ion collisions using one meson 
exchange model with the effective Lagrangian. To take 
into account the effects due to two different spin possi- 
bilities of the X(3872) meson, we have evaluated two ab- 
sorption cross sections for both .P = 1+ and 2~ states of 
the X(3872) meson by pions and rho mesons during the 
hadronic stage of heavy ion collisions. We have found 
that the absorption cross sections and their thermal av- 
erages are strongly dependent on the structure and the 
quantum number of the X(3872) meson; the energy de- 
pendence of the cross sections is quite different for two 
spin states as shown in Fig. [2l and the cross sections 
arc much bigger for a 2~ state than those for a 1+ state. 
Therefore, it is expected that the spin-2 X(3872) meson 
can be absorbed by light mesons much more easily than 
the spin-I X(3872) meson. 

We have further investigated the time evolution of the 
abundances for two possible quantum number states of 
the X(3872) meson. We have found that the variation 
of the X(3872) meson abundance during the expansion 
of the hadronic matter is also strongly affected by the 
quantum number of the X(3872) meson; the Xi(3872) 
meson abundance slightly changes to 4.3x10"^, while 
the X2(3872) meson abundance varies significantly to 
1.2x10"'*, leading to the final abundance ratio ^ 2.8 
between X2(3872) and Xi(3872) mesons at the kinetic 
freeze-out. We therefore suggest that studying the abun- 
dance of the X(3872) meson in relativistic heavy ion col- 
lisions provides a chance to infer its quantum number as 
well as its structure. 

We have also considered the hadronic molecular state 
D°D*^{D^D*°) possibly produced from D^{m) and 
£)*Oj-£)*0'j sometime during the hadronic stage. The 
abundance is expected to be in the range between 
7.8xl0~* and 8.1 xlO""*, resulting in the final ratio ~ 18 
between the hadronic molecular state and the tetraquark 
state. The ratio is large enough to be used to discrimi- 
nate the structure of the spin-1 X(3872) meson. 
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Appendix A: Dependence of the strong coupling 
constants on the X(3872) meson mass 



with 



In order to understand the origin of the big difference 
between gxiD* d and gx2D* d we investigate the relation 
between two strong couphng constants. We simply con- 
sider the decay rate of the X(3872) meson decaying to 
D° and L>*o, 




3871.8 3872.4 3873.0 3873.6 3874.2 



Mx,3S72, (MeV) 



FIG. 6: The strong coupling constants QxiD'd and .gxgiJ'D 
as functions of the X(3872) meson mass. 



f)*0\ 



dT{X D'^D 



1 



(2^)' 



2sx + 1 2mx 



^|Alpd$2(fc;p,g), 



pol 



(Al) 

where (i$2(k;p,q) is the element of 2-body phase space 
density given by 



(A2) 



In the rest frame of the X(3872) meson, the decay rate 
Eq. (|A1[) becomes, 



r(X^D»D-).i--l-H^i:iMP (A3) 



" 2mx 

(A4) 

Since we consider two spin possibilities for the X(3872) 
meson to explain one experimental observation, we ob- 
tain the following condition using the interaction La- 
grangians, Eq. (|H), 



pol 



(A5) 



pol 



by requiring = ^Xi- This condition is responsi- 
ble for the difference between gxiD'D = 3.5 GeV and 
gx2D*D = 189 when a mass of the X(3872) meson is 
3872.26 MeV. We can obtain, using Eq. ^M^, the other 
strong coupling constant of the X(3872) meson when one 
of them is known. 

We have also found that gx^D'D is sensitive to the 
variation of the X(3872) meson mass when obtaining it 
from the requirement Tx^ = Tjfi ■ Since we know that 
a recent measurement of the X(3872) meson mass from 
its decay mode, X^ D*°D^ is 3872.9 MeV [U, we vary 
the mass of the X(3872) meson, and see how the strong 
coupling constants gxiD*D and gx2D*D change. 

As shown in Fig. [SI gx2D*D decreases from 189 to 
107 whereas gxiD*D changes from 3.5 GeV to 2.5 GcV 
when the mass of the X(3872) meson increases from 
3872.3 MeV to 3872.9 MeV. We expect the above new 
strong coupling constants obtaiend when the mass of 
the X(3872) meson is 3872.9 MeV to reduce the cross 
sections evaluated in Sec. Ill by a factor of 3 for the 
X2(3872) meson, and by a factor of 2 for the Xi(3872) 
meson, respectively. 
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